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Abstract 

Water availability is a major environmental driver affecting riparian and wetland vegetation. The interaction between water 
table fluctuations and vegetation in a stochastic environment contributes to the complexity of the dynamics of these ecosystems. 
We investigate the possible emergence of spatial patterns induced by spatio-temporal stochastic resonance in a simple model of 
^-H groundwater-dependent ecosystems. These spatio-temporal dynamics are driven by the combined efl'ect of three components: (i) an 
additive white Gaussian noise, accounting for external random disturbances such as fires or fluctuations in rain water availability, 
(ii) a weak periodic modulation in time, describing hydrological drivers such as seasonal fluctuations of water table depth, and (iii) 
^>-^a spatial coupling term, which takes into account the ability of vegetation to spread and colonize other parts of the landscape. A 
^ suitable cooperation between these three terms is able to give rise to ordered structures which show spatial and temporal coherence, 
^> and are statistically steady in time. 
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p I 1. Introduction 

O Random fluctuations can play an important role in deter- 
O mining the composition and structure of ecosystems |[3] [TSl . 
' Fire occurrences ["^l, fluctuations in rain water availability 
, ^1 ET] [TTl [LLI, soil heterogeneity |33|, and temperature oscilla- 
tions 1321 , are some of the examples of fluctuating environ- 
mental drivers. Evaluating ecosystem response to environmen- 
tal fluctuations is crucial, among other things, to the manage- 
ment and restoration of plant communities (HI [521, the under- 
T:::;j- standing of their ability to sustain biodiversity ||26l|23l[T0l, and 
to predict causes of possible shifts in vegetation composition 

EillSlEIl. 

The efl'ects of noise on dynamical systems has been studied 
CN in a number of stochastic models proposed in recent years |42, 
[141. Noise is typically associated with a "destroying efl'ect", 
^ as a driver of disorganized fluctuations around the stable states 
of the underlying deterministic dynamics. However, a number 
%N of studies have demonstrated the existence of a "constructive 
efl'ect" of noise: ordered spatial structures (i.e. patterns) can 
emerge in some spatio-temporal dynamical systems as a result 
of noisy fluctuations (see, for example, (361 [43l |50l [29l [45l ). 
An increase in noise intensity can produce in these systems a 
counterintuitive, more regular and organized, behaviour both in 
time and in space. 

For example, the cooperation of a random external driver 
with an external periodic forcing is able to induce ordered tran- 
sitions between difl'erent states of a dynamical system. Known 
as stochastic resonance (see, among others, |[T3]|20l[5Tl), this 
phenomenon may occur when a bistable system is disturbed by 
an external random forcing in the presence of a weak periodic 
fluctuation in time. Bistable system exhibit a bistable potential 



with two minima separated by a potential barrier. The periodic 
forcing causes fluctuations in the height of the potential bar- 
rier, without ever causing a transition between the two states. 
It is the external random forcing - if strong enough - that is 
able to drive the system from a stable state to the other; these 
transitions are more likely to occur when the height of the po- 
tential barrier is lowest, which depends on the phase of the pe- 
riodic fluctuations. A suitable synchronization between the fre- 
quency of the random transitions (which is related to the noise 
intensity) and the frequency of the periodic forcing induces a 
resonance-like efl'ect that generates regular transitions between 
the two stable states (53l[3n . 

Even though a number of studies have recognized the 
bistable character of some ecosystems (39| [5] [56l |37l and the 
ability of random fluctuations to induce new dynamical be- 
haviours that do not exist in the underlying deterministic dy- 
namics (21 in, applications of stochastic resonance theories to 
eco-hydrology have started to appear only recently (see, for in- 
stance, |4l|46l[3ll47J). 

In this paper we study the basic mechanisms for the oc- 
currence of spatio-temporal stochastic resonance in a simple 
model describing the dynamics of riparian or wetland vege- 
tation dynamics. The temporal deterministic dynamics of the 
model have been proposed by Ridolfi et al. (37l to investi- 
gate the emergence of bistability from the interactions between 
phreatophyte vegetation and shallow phreatic aquifers in wet- 
land ecosystems. The combined eff'ect of external noise and pe- 
riodic oscillations on these bistable dynamics has been shown 
to be able to induce stochastic resonance in the time domain 

HI 

We investigate the possible emergence of spatio-temporal 
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stochastic resonance when a diffusive spatial coupling term is 
introduced in the presence of random and periodic forcings. 
The addition of a spatial diffusive mechanism accounts for the 
vegetation ability to encroach across the landscape. The other 
assumptions are physically based, due to the presence of envi- 
ronmental disturbances (e.g. occurrence of fires, rainfall fluc- 
tuations) acting on the system, and a number of periodic hy- 
drological drivers (e.g., seasonal oscillations of the water ta- 
ble). In particular, a significant connection exists between wa- 
ter table fluctuations and vegetation dynamics. The water level 
variation appears to be an independent and important vegeta- 
tion gradient, as the responses of species to the range of water 
level fluctuations seem to reflect their tolerance to disturbances 
ifTSl . Changes in species composition and distribution, as well 
as in vegetation structure, have been related to various factors, 
mainly the vegetation topographical position with respect to the 
water table, frequency and duration of inundation ll24ll . It has 
been observed that, under the periodic declining groundwater 
availability, vegetation patches tend to show a regular spatial 
behaviour |[T6ll . Moreover, in dune systems, the occurrence of 
regular patterns of heath patches has been associated, among 
other things, to the different groundwater discharge |25|. 

The relevant question here is whether the phenomenon of 
stochastic resonance - resulting from the combined effect of 
noise and a weak temporal oscillation - can interact with a spa- 
tial diffusive term to generate vegetation ordered spatial struc- 
tures. 




Figure 1 : Vegetation carrying capacity Ycc as function of the vegetation biomass 
y, with = 11 and d^up = 1.8. (a) = 1, dinf = 0.95. (b) = 1, dinf = 1.2. In 
the upper panel the potential is shown, (c) ^ - 0.6, dmf - 1.2. 

The depth, d, of the water table is adimensionalized with re- 
spect to water table depth in the absence of vegetation, J*, and 
can be expressed as a linear function of plant biomass, V, 

J=l+ySy, (2) 

where is the dimensionless sensitivity of the water table to 
the presence of vegetation, is positive because plants typi- 
cally tend to increase the depth of the water table. The water 
table depth, J, in turn affects the dynamics of wetland vegeta- 
tion. If the water table is too shallow (waterlogging), vegetation 
can suffer, due to an insufficient aeration of the root zone and 
a decreased rate of seedling establishment BTl . If the water 
table is too deep, water is out of reach of taproots and vegeta- 
tion can suffer as well. These effects are taken into account by a 
quadratic dependence of the carrying capacity, V^c, on the water 
table depth, J, 



2. Modeling framework 

We consider the dynamics of plant biomass, V, and its inter- 
play with the phreatic aquifer |37]. We first recall some impor- 
tant aspects of the temporal deterministic dynamics, and then 
present the characteristics of the spatio-temporal stochastic sys- 
tem. 

2.7. Deterministic model 

Changes in vegetation biomass are the result of a growth- 
death process that can be expressed as (e.g., 128. 48J ), 

dV 

= V(Vcc - V), (1) 

dr 

where V is the dimensionless biomass, normalized by a (high) 
fixed reference value, r = at is the dimensionless time (where a 
controls the temporal response of the system and t is time), and 
Vcc is the dimensionless ecosystem carrying capacity, that is the 
maximum amount of vegetation sustainable with the available 
resources. 

The interaction between phreatophyte vegetation, i.e. plants 
relying on the phreatic aquifer, and the average depth of the lo- 
cal water table is widely recognized as one of the key aspects 
affecting wetland ecosystems dynamics (301 SD H [37l . For ex- 
ample, the water table decreases in the presence of phreatophyte 
species because of the lower recharge rates due to rainfall inter- 
ception and plant transpiration |[54j[T]|38l[T2l, and also because 
vegetation taproots directly extract water from the aquifer |[T9l . 



a(d - dinf){dsup - d) if dinf <d < 
1 otherwise, 

where a regulates the sensitivity of carrying capacity to changes 
in water table depth, dinf and dsup are the adimensional thresh- 
olds of vegetation tolerance to shallow and deep water tables, 
respectively. 

A positive feedback may exist between vegetation establish- 
ment and water table dynamics, whereby phreatophytes favor 
their own survival by increasing the water table depth and en- 
hancing root aeration [54, 7]. The presence of positive feed- 
back mechanisms suggests the possible existence of multiple 
stable states |37|. This fact is confirmed by some experimen- 
tal evidence suggesting that two alternative equilibrium states 
may exist in wetland vegetation dynamics ||4TJ|6l[55l|44l. Equi- 
librium states can be found by setting Eq. ([T]) equal to zero, 
y = Vo = is always an equilibrium state, therefore the exis- 
tence of multiple equilibrium states depends on the existence of 
real, nonnull roots of the equation Vcc = V- These solutions can 
be found by substituting Eq. ^ in Eq. ([3]), and then search- 
ing for intersections between VcdV) and Vcc = V (see Fig. [T]). 
Ridolfi et al. [37 1 showed three possible cases: 
(a) dinf < 1 (Fig. [T^). The water table depth in the absence of 
vegetation, d = 1, is always greater than the minimum depth 
required for vegetation establishment. Vq = is an unstable 
state and only one stable state exists, Vi . Ecosystem dynamics 
always tend the vegetated state, Vi , regardless the initial condi- 
tions. 
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(b) dinf > 1 (Fig. [T|)) and vegetation is able to keep the water 
table below the maximum elevation acceptable for plants (e.g. 
P = I), three equilibrium states are found. Two of these states, 
y = Vo and V = Vi are stable, while the other one, V = V^. is 
unstable. If V grows above V^, then the system tends to the veg- 
etated state. Otherwise, it evolves towards the unvegetated state 
and remains here blocked. Note that the two stable states cor- 
respond to the minima of the double- well potential, U(V), with 
^ = - ^ , while the unstable solution, V = corresponds to 
a maximum or "potential barrier". 

(c) dinf > 1 (Fig. [T]:) and vegetation is not able to reduce the 
water table elevation below the minimum level suitable for plant 
survival (e.g. p = 0.6). In this case the only possible stable state 
is the unvegetated condition, V = Vq = 0. 

2.2. Spatio-temporal stochastic model 

We now focus on the case (b) described above, with two sta- 
ble states (vegetated and unvegetated) emerging from the de- 
terministic dynamics. We introduce in these dynamics three 
new terms: (i) a periodic temporally oscillating term, account- 
ing for seasonal or interannual fluctuations of the water table 
depth, J, 1401 [221 ; (ii) a random forcing term, representing the 
efl'ect of random environmental drivers (e.g., fires or climate 
fluctuations); (iii) a spatial coupling term, modeling the diflfu- 
sive spread of vegetation. 
Eq. ([T]) now reads. 



dv 



= v(y,,-y) + ar,T) + z)v2y 



(4) 



where D is the strength of the spatial coupling, is the Laplace 
operator, ^(r, r) is a white (in time and space) Gaussian noise 
with zero mean, intensity s, and correlation (^(r, T)^(r',TO) = 
2sS(r - r')S(T - r'), being S(') the Dirac delta function. The 
temporal oscillation is directly inserted into Eq. (|4L through 
the dependence of Vcc on the water table depth, EqTQ, which 
can be now expressed as 



d = 1 -\- Acos 



(5) 



where A and oj = y the nondimensional amplitude and 
the dimensional frequency of the seasonal oscillations, respec- 
tively. T is the oscillation period (for example T = I year). No- 
tice that these seasonal fluctuations of the water table are crucial 
to the emergence of the dynamical behaviors (i.e., spatiotempo- 
ral stochastic resonance) presented and discussed in the follow- 
ing sections. Therefore, Eq. ^ now reads 



[l + Acos (f ) +/5y - dinf] [dsup - 1 - Acos (f ) -fiv] 







if dinf <d <dsi 
otherwise. 



(6) 

We consider the case of relatively weak seasonal oscillations 
in a way that the potential, UiV), still has two minima corre- 
sponding to two stable states. Fig. [2]shows the potential, UiV), 
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Figure 2: Potential U{V) of the temporal deterministic model ^ at different 
times during the periodic cycle of d jsj. A = 0.08, a = 0.5 d,co = 0.0344 
rad/d, a = 13, dinf = 1.2, d^up = 1.8,/? = 1. 

of the temporal deterministic model ([T]), but with the water ta- 
ble depth, J, modulated as in ([5]). The potential always main- 
tains two minima throughout its periodic oscillation. Thus, even 
though there are phases in the modulation of the periodic forc- 
ing in which the height of the potential barrier is strongly re- 
duced (hence, the probability of occurrence of a state transition 
is enhanced), deterministic transitions between the two wells of 
the potential are not allowed. 

The choice of an additive white (Gaussian) noise to repre- 
sent external drivers acting on vegetation dynamics in the time 
domain is motivated by the fact that the temporal scales of the 
random driver are typically much shorter than the characteristic 
temporal scales over which the temporal vegetation dynamics 
evolve (e.g. the relaxation time to a stable state). Moreover 
we used a white noise also in space to show how spatial co- 
herence is not imposed by the spatial correlation of these exter- 
nal driver but emerges as an efl'ect of spatio-temporal stochastic 
resonance. Therefore, this kind of noise is typically adopted in 
stochastic modeling [i42il . 

The Laplacian, V^, in Eq. Q is a simple operator which 
is widely used to represent the spatial eff'ects of the difl'usion 
mechanisms in vegetation dynamics |3l|22|[52l. This operator 
accounts for spatial interactions between a point of the domain 
and its nearest neighbors, and is therefore considered as a short- 
range spatial coupling. 

3. Results 

The model ^ is numerically solved for difl'erent parameter 
sets to directly investigate the associated spatio-temporal veg- 
etation dynamics. Periodic boundary conditions are adopted, 
while, in the absence of other indications, the initial conditions 
used in the simulations are uniformly distributed random num- 
bers in the range [0.29, 0.31]. We also use the mean field anal- 
ysis, a mathematical framework that provides the steady- state 
probability distribution of y |[36l |49l [50l . Details about the 
numerical scheme and the mean-field analysis are reported in 



[Appendix A 

Fig. |3] shows numerical results for s = 0.008 and D = O.l 
(other parameter values are reported in the caption of the fig- 
ure). Black and white tones are used for the highest (V = 0.7) 
and lowest (y = 0) values of V, respectively. Well defined 
patterns periodically emerge (e.g., at r = 250 and r = 341) 
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Figure 3: Model s = 0.008, D = 0.1, A = 0.08, a = 0.5 d'^ = I, 



13, d,, 



1.8, d, 



inf 



1.2, T = 182.5 days. Top: numerical simulations 



at T = 250, 295, 341. Bottom: mean and variance of the vegetation biomass, V, 
as functions of time (black circles correspond to mean and variance values at 
T = 250, 295, 341). Pdf of vegetation field (solid: mean-field analysis, dashed: 
numerical evaluation). 



and disappear (r = 295) in the course of periodic fluctuations 
of d. In this case the non-dimensional period is T' = 91.25. 
Pattern occurrence corresponds to temporal minima of mean 
vegetation biomass and to maxima of biomass variance, while 
homogeneous vegetated states correspond to maxima of mean 
biomass and minima of the variance of vegetation biomass (see 
circles in Fig. [3]). Here and in the following discussion, we call 
"homogeneous" a state in which the system exhibits negligible 
spatial heterogeneity with respect to the patterned state. 

The spatio-temporal dynamics continuously oscillate be- 
tween two difl'erently vegetated states. In this example patterns 
appear during the transition between these two states, when the 
system approaches the least vegetated stable state, and disap- 
pear when the most vegetated state is reached. The interaction 
of noise with periodic water table oscillations and spatial dif- 
fusion gives rise to spatial ordered structures which are statisti- 
cally steady and periodically appear/disappear in time with the 
same period, T, as the fluctuations in d. These periodic state 
transitions, and the associated formation and disappearance of 
spatial patterns are the evidence of spatio-temporal stochastic 
resonance, with noise-induced shifts between the two minima 
of the potential, U(V). Although the system tends to oscillate 
between two stable states, these states are not always reached. 
Indeed, temporal oscillations have to be sufficiently long to per- 
mit the system to visit both states. Here, for example, the tem- 
poral period, T = 182.5 days, is too short and does not allow 
the system to reach the unvegetated state. Conversely, with a 
longer period, T, the system is able to attain both the vegetated 
and the unvegetated states (see Section [33] ). 

In Fig. [3] numerical and analytical evaluations of the pdf 
of the vegetation field are reported. Here and in the follow- 
ing discussion, the numerical approximation of the pdf (dashed 
curves) is obtained using simulated fields in the temporal range 
50 < T < 547, which encompasses several temporal oscillation 
periods (here, T' = 91.25). Details on the analytical approxi- 
mation of the pdf (solid curves) are discussed in [Appendix A 





The agreement between numerical and analytical evaluations of 
the pdf is quite good in showing a weak bimodality. 



Figure 4: Numerical simulations of model (|4j with D = 0.2, other parameters 
as in Fig. |3] (Top row) Field distributions with s = 0.0024 at r = 275 (left 
column), T = 320 (middle column), and pdf (right column). (Middle row) Field 
distributions with s = 0.012 at r = 202 (left column), r = 247 (middle column), 
and pdf (right column). (Bottom row) Field distributions with s = 0.04 at 
T = 270 (left column), r = 310 (middle column), and pdf (right column). 



3.1. Role ofD and s 

We now consider changes in the spatial coupling strength, D 
(D = 0.01, 0.1, 1), with a fixed noise intensity (e.g. s = 0.008). 
As D increases, the least vegetated state moves closer the most 
vegetated one. As a consequence, for higher D values, the sys- 
tem becomes more homogeneous (the biomass variance oscil- 
lates around lower values and the amplitude of these fluctua- 
tions becomes smaller). By increasing D from 0.01 to 0.1, spa- 
tially coherent patterns appear. If D is further increased, the dif- 
fusive efl'ect is so strong that the field is almost homogeneous, 
and the system experiences a blocking eff'ect induced by diff'u- 
sion. 

Diff'erent noise levels (s = 0.0024,0.012,0.04) with a fixed 
spatial coupling (D = 0.2) are then considered. If the noise 
intensity is weak (s = 0.0024) with respect to the spatial cou- 
pling strength, D, the biomass density tends to be homogeneous 
in time and space. On the other hand, if the noise level is too 
strong (s = 0.04), the system experiences random oscillations 
between more vegetated and less vegetated states. However, in 
this case noise has a "destroying" eff'ect in that it destroys any 
spatial and temporal coherence and no ordered states emerge. 
An intermediate noise value (s = 0.012) allows for the emer- 
gence of temporal fluctuations and the formation of periodic 
spatial patterns. The existence of an intermediate noise level 
which, in cooperation with a weak temporal oscillation, is able 
to regularly move the system between difl'erent states is the 
most typical feature of the stochastic resonance phenomenon 
(53] [131. In Fig. [4] we show the configuration of the system 
in correspondence to local minima (r = 275 for s = 0.0024, 
T = 202 for s = 0.012, r = 270 for s = 0.04) and maxima 
(r = 320 for s = 0.0024, r = 247 for s = 0.012, r = 310 
for s = 0.04) of biomass variance (left and central columns. 
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Figure 5: Model s = 0.0024, D = 0.2, other parameters as in Fig. |3] (Top 
row): Field distributions with random initial conditions uniformly distributed 
in the range [0, 0.02] at r = (left column), r = 250 (middle column), and 
pdf (right column). (Bottom row): Field distributions with random initial con- 
ditions uniformly distributed in the range [0.29, 0.31] at r = (left column), 
T = 250 (middle column), and pdf (right column). 



Figure 6: Model jJ), s = 0.008, D = 0.1, A = 0.08, a = 0.5 6-^/3= 1, 
a = 13, dsup = 1.8, dinf = 1.2, T = 365 days. Numerical simulations at 
T = 235, 277, 310, 342. Bottom: mean and variance of the vegetation biomass, 
V, as functions of time (black circles correspond to mean and variance values 
at T = 235,277, 310, 342). Pdf of vegetation field (solid: mean-field analysis, 
dashed: numerical evaluation). 



respectively) for the three values of noise intensity. 

In general, temporal local maxima (minima) of biomass vari- 
ance correspond to minima (maxima) of biomass mean only 
when the system is able to exhibit ordered spatial structures 
(see Fig. |3]). In all the other cases - i.e. with homogeneous 
(Fig. |4] top) or disturbed (Fig. [4j bottom) spatial distributions - 
there is not such a synchronization between maxima (minima) 
of biomass variance and minima (maxima) of biomass mean. 

Also in these cases, there is a very good agreement between 
the mean-field steady-state pdfs (solid curves) and the pdfs nu- 
merically calculated (dashed curves) using simulated fields in 
the temporal range 50 < r < 547. This result is especially in- 
teresting because this mean-field analysis accounts for the mod- 
ulation of the dynamics in time (see [Appendix A ), an approach 
that has never been reported before. 

From the results here presented, we can conclude that there 
is an optimal region in the (s, D) plane, in which spatial patterns 
emerge. This region corresponds to a noise intensity, of order 
10"^, and a spatial coupling strength, Z), of order of magnitude 
10-1. 

3.2. Influence of initial conditions 

In the previous sections the initial conditions were generated 
as random numbers uniformly distributed in the interval [0.29, 
0.31]. If noise is sufficiently strong with respect to the spatial 
coupling strength, Z), results similar to those shown in figures 
[3] and |4] (middle and bottom), are obtained with diff'erent initial 
conditions. Conversely, if the noise intensity, s, is small with 
respect to Z), the system tends to an homogeneous steady state, 
which depends on the initial conditions. In other words, a suffi- 
ciently high noise intensity with respect to the spatial coupling 
strength allows the system to continuously oscillate in time be- 
tween the two minima of the potential, regardless the initial 
conditions. On the contrary, if noise is small with respect to 
the spatial coupling strength, the system remains trapped into 
the minimum of the potential which is closer to the initial con- 
ditions. 



An example is shown in Fig. [5] where, for s = 0.0024 and 
D = 0.2, diff'erent initial conditions, uniformly distributed in the 
interval [0, 0.02] or [0.29, 0.31], lead to an unvegetated (top) or 
a vegetated (bottom) state, respectively. When the initial con- 
ditions are in the interval [0, 0.02], the system is initially close 
to the minimum of the potential (Vq = 0). The noise intensity s 
is so weak with respect to the spatial coupling strength, D, that 
it is not able to move the system away from this stable state; 
therefore the system remains blocked into the unvegetated sta- 
ble state. 

3.3. Difl'erent periodicity of the temporal modulation 

We then evaluate the eff'ect of a diff'erent temporal periodic- 
ity in the modulation term, G = cos(^), by taking T = 365 
days, while all the other parameters remain the same as in Fig. 
[3] Patterns periodically emerge (r = 277 and r = 342) and dis- 
appear (r = 235 and r = 310). Their occurrence corresponds to 
temporal local maxima of the vegetation variance, while homo- 
geneous vegetated and unvegetated stable states correspond to 
maxima and minima of the vegetation mean, respectively (see 
circles in Fig. |6]). A longer period, T, allows the system to ex- 
plore the whole range of values of V. This fact is confirmed 
by a non-negligible bimodality of both analytical and numeri- 
cal pdfs, whose peaks correspond to the two alternatively vis- 
ited stable states. During these state transitions, patterns appear. 
Although the temporal modulation establishes the frequency of 
pattern appearance, the mechanism of pattern formation is al- 
ways present and does not depend on the specific temporal pe- 
riodicity imposed into the system. 

3.4. N on- simultaneous presence of random and periodic forc- 
ings 

We now consider the case in which noise and temporal mod- 
ulation do not occur at the same time. As noted, as the noise 
level, s, decreases the spatial distribution of plant biomass be- 
comes more homogeneous. A similar behavior is expected to 
occur for ^ = 0. The homogeneous state depends on the initial 
conditions: if the system is initially in the interval [0, 0.02] it 
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Figure 7: Model (|4| with s = 0.008, A = 0, other parameters as in Fig. |3] 
initial conditions are random numbers uniformly distributed in the interval [0, 
0.02]. (Top row) Field distributions with Z) = 0.1 at r = 0, 125, 250 (from left 
to right), and pdf (right column, dashed: numerical at r = 547, solid: mean- 
field analysis). (Bottom row) Field distributions with D = 1 at r = 0, 125,250 
(from left to right), and pdf (right column, dashed: numerical at r = 547, solid: 
mean-field analysis). 



homogeneously decays to zero, while for initial conditions in 
the interval [0.29, 0.31] it reaches the homogeneous vegetated 
stable state and weakly oscillates in time around it. 
When noise is present with no modulation (A = 0), state transi- 
tions may occur. If the noise is weak with respect to the strength 
of the spatial coupling, D, the homogeneous steady state at- 
tained by the system depends on the initial conditions. If the 
noise is strong enough with respect to D, the system tends to 
the homogeneous vegetated state regardless of the initial condi- 
tions. Patterns can emerge only as transient features before the 
system approaches the homogeneous steady state. An example 
is shown in Fig. [7] where ^ = 0.008, the initial conditions are in 
the interval [0, 0.02], D = OA (top) and Z) = 1 (bottom). For 
D = 0.1, while passing from an unvegetated initial condition 
(r = 0) to a homogeneous vegetated steady state (r = 250), the 
system exhibits some transient patterns (r = 125), which fade 
out in the long run. Instead, for Z) = 1 the noise is weak with 
respect to D and therefore the system is blocked into the low 
vegetated state without even showing transient patterns. 



4. Discussion and conclusions 

In the system investigated in this paper the occurrence of 
spatio-temporal stochastic resonance can be thought of as the 
result of the combined effect of three components: a white (in 
time and space) Gaussian noise, a weak periodic modulation 
in time, and a suitable spatial coupling term. In the absence 
of modulation, the noise - if strong enough with respect to the 
spatial coupling - can lead the system towards the vegetated 
state, regardless of the initial conditions. However, if the noise 
is weak with respect to D, the system will approach one of the 
two (vegetated or unvegetated) stable states, depending on the 
initial conditions. With no modulation, patterns can only occur 
in an initial transient because, once the stable state is reached, 
the system is locked in a homogeneous state. On the other hand, 
in the absence of noise the temporal modulation is not able to 
promote a state transition, because the system reaches one of 
the two stable states - depending on initial conditions - and 
weakly oscillates around it. The spatial coupling term is able 
to induce spatial coherence without itself enhancing state tran- 
sitions. In other words, noise is the driver able to induce state 



transitions, while the periodic forcing and the spatial coupling 
are able to induce temporal and spatial coherence, respectively. 
Only a suitable cooperation between these three mechanisms 
can lead to the interesting ordered scenarios analyzed in Sec- 
tion [3] Thus, random and periodic forcings have to be simulta- 
neously present along with a suitable spatial coupling to induce 
the formation of statistically steady ordered structures with spa- 
tial and temporal coherence. 

Changes in the periodicity, T, of the weak temporal modu- 
lation strongly affect vegetation dynamics. In this case, if T is 
of the order of six months or one year, coherent spatial struc- 
tures periodically appear during transitions between vegetated 
and unvegetated states. These dynamics could represent the 
emergence and disappearance of an herbaceous species, which 
typically have a growth/life period comparable to the driving 
seasonal oscillations of the water table. Conversely, the emer- 
gence of patterns in model simulations with longer periods may 
represent the formation of ordered structures in the distribution 
of tree species, which grow over longer time scales. 
Beside this interpretation based on a single species, spatial pat- 
terns result from the simultaneous existence of two different 
species, whose growth/life periods are comparable to T but not 
synchronized with one anther (i.e. when species 1 is at max- 
imum density, species 2 is at minimum density). The first 
species is dominant at a certain time (almost homogeneous 
field, see, for example, left top panel of Fig. [6]), but it tends to 
disappear, while the density of other increases. The simultane- 
ous presence of both may give rise to ordered spatial structures 
(spatial pattern occurrence, see middle top panel of Fig. [6]). At 
some point, the second species temporarily dominates the plant 
community composition (almost homogeneous field, see right 
top panel of Fig. [6]), but it will then start to decay, while the 
first species grows back and new patterns emerge (see left bot- 
tom panel of Fig. |6]). 

Appendix A. Numerical and analytical methods 

The typical numerical approach used to solve stochastic par- 
tial differential equations is based on a discretization of the con- 
tinuous spatial domain using a regular Cartesian lattice with 
spacing Ax = Ay = A. A two-dimensional square lattice with 
128x128 sites and A = 1 is here adopted. The original equa- 
tion ^ is then transformed into a system of coupled stochastic 
ordinary differential equations with a finite difference scheme. 



dVi 
dr 



= ViiVcci - Vi) + ^i + D (Vj - Vi)^ (A. 1) 



jenn(i) 



where Vi and ^/ are the values of V and ^ at site /, respectively, 
nn(i) is the set of the four nearest neighbors of the site /. Nu- 
merical simulations are carried out with the Heun's predictor- 
corrector scheme |[50l|42l, with a temporal step Ar = 5 • 10"^. 
Notice how this method does not prevent the emergence of a 
spurious spatial correlation associated with the discrete repre- 
sentation of the dynamics in a 2D lattice | 21 ]. In our analyses 
the simulations were repeated with grids of different sizes and 
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the same qualitative patterns and dynamical behaviors were ob- 
served to emerge regardless of the resolution used in the dis- 
crete representation of the domain. 

A qualitative representation of (stochastic) spatio-temporal 
dynamics can be obtained through the mean-field method. Its 
fundamental assumption is that the mean of the values of V in 
all the neighboring cells can be approximated by the spatio- 
temporal mean of the field, namely Z/^;) = 4(Vy) = 4(y). 
Therefore, eq. (|A.1|) becomes 



dr 



Vi(Vcc,i-Vi)^^i^4D({V}-Vi). 



(A.2) 



The classic mean-field analysis allows one to evaluate if a 
phase transition occurs. The existence of multiple solutions, 
{V}n, of the self-consistency equation 



{V}n 



f 

Jo 



Vp'J(V\G;{V)J dV = F{{V)nl 



(A3) 



indicates the presence of a non-equilibrium phase transition. 
In Eq. {A3) G = cos(^) is the temporal modulation and 
p^J(V\G; {V}n) is the conditional probability with respect to G at 
steady state. We recall that the occurrence of non-equilibrium 
phase transition is neither a necessary nor a sufficient condi- 
tion for noise-induced pattern formation |[36l|43l. In fact, non- 
equilibrium phase transitions imply that noise is able to change 
the value of the order parameter, but not that ordered geometri- 
cal structures necessarily emerge. 

Here, the dynamics exhibit also a temporal modulation, G, 
which has to be adequately accounted for while applying the 
mean-field method. To evaluate the steady state pdf of the 
biomass vegetation, V, and at the same time to deal with 
the temporal modulation G, we propose to solve the self- 
consistency condition ( |A.3| ) for all the values of G e [-1,1]. 
In this way the conditional probability, p^^XV\G; <y)„), is ob- 
tained. The steady-state probability p^XV', {V}) is then com- 
puted through the convolution product between the conditional 
probability, Pn(V\G; <y)„), and the probability, p(G), of the 
temporal modulation, G. 

Depending on s, D and G values, one or three solutions of 
the self-consistency equation (|A.3|) can be found. If the solu- 



tion of ( |A.3| ) is unique, then it is considered as the steady state 
solution. If three solutions are obtained, we can assume that the 
lowest one corresponds to the least vegetated stable state (Vq), 
the highest one to the most vegetated stable state (Vi), while 
the one in the middle to the unstable state (Vu). If the temporal 
period, T, is sufl&ciently long (e.g. T = 365 or 730 days) to al- 
low the dynamics to fully develop and periodically visit the two 
stable states, all the three solutions are taken into account as 
steady state solutions. If, instead, the period T is much shorter 
(e.g. T = 182.5 days), the dynamics are able to completely 
reach only one of the two stable states, while the other one re- 
mains unexplored. As a result, the system oscillates in time 
and fully captures only one stable state, while the other one is 
not approached. Only the explored stable state and the unstable 
state are considered as steady state solutions. 



Once solutions, {V}n = of ( |A.3| ) are determined, their 
conditional probabilities, p^J(V\G', {V}n), are computed as 



p'JiVlG; <y)J = - cxp[(-U(V) - 2DV^ + 4D/j,V)/s] (A.4) 
and 



Z, = r exp[(-f/(y)-2W" + 4D//,y)/^]dy, (A.5) 
Jo 

where f/(y) is the potential of the temporal deterministic model 
([T]) with water table depth, J, modulated as in ([5|. Since the pdf, 
p(G), of the cosine function, G = cos (^), is 



piG) 



1 



(A.6) 



the probability p^^{y\ (y)) can be evaluated through the convo- 
lution integral 



/'(V;<y)) 



(V|G;(y)J 



p(G)6G. (A.7) 



If, for a given value of G, two or more conditional probabili- 
ties, p^^{y\G\ (y)«), exist, they are summed together and then 
weighted with the same probability (defined by p(G)) to com- 
pute the probability p'^iy\ (y)). The probability p'^iy\ (y)) is 
in the end normalized so that p^^{y\ <y))dy = 1. 
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